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Abstract. We analyze ISO/SWS spectra of the RV Tau star R Scuti. The infrared spectra are dominated by 
H2O emission bands. The near- and mid-infrared excess is attributed to H2O; the dust contribution is less 
important. We also identify CO, SiO and CO2 bands. The various molecular emission bands originate from an 
extended atmosphere, an atmosphere above the photosphere. The extended atmosphere of R Set is formed from 
matter which gradually have lifted up from the photosphere through the pulsations of the star. In contrast to 
the abundant molecules around the star, the silicate dust feature is weak and the dust mass-loss rate is only 
Afd = 10 -11 MQyr -1 . This implies that there might be a process to inhibit dust formation from molecules. RV 
Tau stars are commonly considered as post-AGB stars. While a detached dust envelope around R Set is consistent 
with such an interpretation, we show that its period evolution is slower than expected. We argue that R Set may 
be a thermal-pulsing AGB star, observed in a helium-burning phase. 
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1. Introduction 

RV Tau stars are pulsating variables characterized by al- 
ternating deep and shallow minima in their light curves. 
They are generally considered to be post-AGB stars with 



found in R Set from the analysis of sodium lines (Lebre 
& Gillct 1991 ). High- resolution spectroscopic observations 
have uncovered CO, H 2 Q a nd OH molecules in the K-band 
(Mozurkewich et al. 1987). The complicated variation of 



low initial masses (~1Mq; Jura 1986). The abundance 

I 



ratios show that they have experienced first dredge-up at 
the bottom of the red giant branch (Giridhar et al. 2000| ) . 
Based on their infrared dust excesses, RV Tau stars are 
classified into two groups: those with extensive warm dust 
and those without evidence of dust in the near-infrared 



region (Goldsmith et al. |1987| ). R Set is the brightest star 
in the visible in the latter group. 

R Set (HR 7066; HIP 9220 2) has a reported period of 
147 days (Kholopov et al. 1988 ). The effective t emper ature 
varies from 4750 to 5250 K (Shenton et al. |1994|); the 



the velocity profile shows that the CO originates in the 
upper atmosphere, above the photosphere, which is af- 
fected by the pulsation. Since the spectral-type of R Set 
is too early for the presence of H2O, they conclude that 
the observed H2O is circumstellar origin. 

In this paper, we analyze spectra of R Set taken with 
the Shor t- Wavelength Spectrometer (SWS; de Graauw et 
al. 1996 ) on board the Infrared Space Observatory (ISO; 



spectral type may vary as late as M3 at minimum phase 
(Kholopov et al. |1988[ ). 

An atmosphere beyond the photosphere, which we call 
an extended atmosphere later on in this paper, has been 
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Kessler et al. 1996 ). Because there is no interference from 
molecules in the terrestrial atmosphere, ISO spectra are 
ideal for the study of molecules in stars. We find that the 
IR spectra of R Set is dominated by molecular emission 
features, especially from H 2 0. In addition, SiO and CO2 
bands are identified. We discuss the distribution of these 
molecules and the mass-loss history. The 200-year light 
curve is used to address the evolution of R Set. 

2. Observational data and analysis 

We obtained the SWS data for R Set from the ISO Data 
Archive. R Set was observed on 10 March, 1996 by pro- 
gramme of MOLBANDS (P.I. A. Heske). The phase is 
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Fig. 1. The visual light curve of R Set (Mattei [1995 , 
private communication) . The magnitudes are averaged in 
5-day bins. The date at which the ISO/SWS observed this 
star is shown by the vertical line. 



estimated as 0=0.60 from the AAVSO light curve (Fig. |l|; 
Mattei 1995, private communication), where = 0.0 cor- 
responds to the nearest previous deep minimum. The 



spectra were obtained using the full-grating scan mode 
(AOT 01, scan speed 1). The wavelength range covers 
2.35-45.2 /im, but the signal-to-noise ratio is poor above 
20 /im and these data are not used. The spectral reso- 
lution is A/AA = 300-500, depending on the wavelength. 
The data were reduced using the SWS Interactive Analysis 
package. The calibration parameters of October 1999 were 
used for the wavelength, detector responsivity, and abso- 
lute flux calibrations. The spectra were re-gridded to a 
constant wavelength resolution of A/AA = 300. 

To identify the molecular features, we calculate molec- 
ular spectra using a circular slab model, which was 
used previously for an alysin g spectra of Mira variables 
(e.g., Yamamura et al. 1999 ). We adopt a multiple layer 
model. This is an approximation for a shell that shows 
a temperature and density gradient; the model layers 
do not necessarily represent discrete shells. Line lists 
for H 2 0, CO2 , SiO and CO were taken from P artridg e 
& Schwenke ( |1997| ), HIT RAN (Rothman et al. |2001a| ), 
Lang hoff fc B auschlichcr ( |1993| ), and HITEMP (Rothman 
et al. 2001b ), respectively. The solar isotopic abun dance 
ratio is adopted for H2O and SiO. Bujarrabal et al. ( 1990 ) 
derive 12 C/ 13 C SilO and we assume 10 % of carbon atoms 
are 13 C. For the energy level populations of the molecules, 
local thermodynamic equilibrium (LTE) is assumed. A 
line width of 5 kms -1 is assumed for all the molecules. 
We use a black body as the background source illuminat- 
ing the molecular layers. We assume that the background 
is a black-body with a temperature of 5000 K. The near- 
and mid-infrared spectrum, which is in the Rayleigh Jeans 
region of the black body, is insensitive to the black-body 
temperature. The synthesized spectrum is not significantly 
affected by the black-body temperature with in the range 
5000-7000 K. 



Table 1. The parameters of synthesized molecular spec- 
tra. We adopt two H 2 components: one is needed to re- 
produce the absorption at 2.7 /im and the other causes 
the strong emission features over the whole wavelength 
region. The size of the water layer, R, is given relative 
to the stellar radius (i?*). The features are indicated as 
either emission (em) or absorption (ab). 
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R 


Appearance 




[K] 


[crn~ 2 ] 


[R*] 




CO 


4000 


3.0 x 10 2i 


1.6 


em 


H20-1 


2200 


1.2 x 10 21 


2.8 


em 


SiO 


= Th 2 o-i 


8.0 x 10 20 


=^?H 2 0-1 


em 


H20-2 


1200 


5.0 x 10 19 


5.5 


ab (2.7 /im) 
em (6.2 /im) 


co 2 


800 


3 x 10 18 


10.5 


ab (4.2 /jm) 
em (13 /im) 



3. Results 

The ISO/SWS spectra of R Set are shown in Fig. |. Below 
4.5 /im, the spectra are similar to those of a typical AGB 
star, Mira (o Cet; spectral type of M5e-M9e), which is 
dominated by water-vapour emission (Yamamura et al. 
1999). Water- vapour bands are mostly seen in emission in 
the near- and mid-IR spectra of R Set, except at ~2.7 /im 
where it is in absorption. CO2 is seen in emission at 13- 
16 /im and in absorption at 4.2 /im. CO at 2.3 /im and 
at 4.6 /im is seen in emission. SiO emission features are 
detected at 4.1 /im. A weak silicate dust excess is present 
around 10 /im. 

The ISO spectra of R Set are compared with the syn- 
thesized spectra in Fig. |^(A) (see also Fig. ||). The spectra 
are well fitted with the parameters listed in Table [ll The 
global shape of the spectra are due to H2O (Fig. y(B)). 
To demonstrate that the water features are actually in 
emission, we also plot separately the background source (a 
5000 K black body). The H2O features can be fitted with 
two components (H2O-I, H2O-2): one component gives 
the emission while the other component contributes the 
absorption at 2.7 /im (but also produces weak emission 
features above 6 /im). 

In Fig. ^, each molecular component of a synthesized 
spectra (CO, SiO and CO2) are compared with the ISO 
spectra. By comparing the observed and synthesized spec- 
tra, the observed features can be identified with H2O, SiO, 
CO, and CO2 bands either in absorption or in emission. 
SiO and CO2 are identified in R Set for the first time. We 
assumed the same excitation temperature of SiO as H2O- 
1 in the synthesized spectra. A weak absorption feature 
near 3.1 /im cannot be explained by H2O (Fig. ||(A)) and 
is also not found in o Cet. This feature can be attributed 
to either HCN or C2H2. HCN is tentatively dete cted in 
the millimeter wavelengths (Bujarrabal et al. 
could be a candidate for the 3.1 /im feature. 



and it 



Radiation from the dust grains was included in the syn- 
thesized spectra. The dust shell was assumed to surround 
the molecular extended atmosphere; the dust has a r~ 2 
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Fig. 2. The ISO/SWS spectra of R Set (black line). The locations of the major molecular features on R Set are 
indicated. The spectra of R Set is compared with that of an oxygen-rich AGB star, o Cet (grey line; flux is scaled). 
These two stars show similar spectra below 4.5 /xm. Note that the spectral resolution of o Cet is twice as high as that 
of R Set. 



density distribution as would be expected for a constant 
mass-loss rate. The temperature of the dust grains was cal- 
culated by solving the energy equilibrium, assuming the 
optically thin condition. For the calculation of the dust 
temperature distribution, the central star was assumed to 
be a black body of 5000 K. The star was then replaced 
by the molecular model spectra, and the radiation trans- 
fer was re-calculated. We used the dust opacity of 'warm 
oxygen-deficient circumstellar dust' (set 1) in Ossenkopf 
et al. fl!992| ). 

The parameters required to reproduce the 10 /im sil- 
icate emission are a dust mass- loss rate of Md = 1-5 x 
10 _11 M^yr _1 , and an inner radius of 50 i?*, assuming 

19881 ) . 



the outflow velocity is lOkms 1 



(Bujarrabal et al. 



The stellar radius (R*) was taken as 6 x 10 cm, and the 
luminosity as 4000 Lq. The outer radius of the dust shell 
in Fig. H| is R a = 5 x 10 4 i?*. However, the spectra can 
be well reproduced within the range of R a — 10 3 -10 5 i?*, 
and we do not well constrain R a . Note that the radii of 
the molecular layers are measured relative to the stellar 
radius; the assumed value of i?*=6 x 10 12 cm is only used 
for the dust model. 

4. Discussion 

4.1. Evidence for molecules in the extended 
atmosphere 

Several molecules are detected in the ISO/SWS spectra of 
R Set. The range of excitation temperatures (800-4000 K) 
is lower than that of the underlying star, confirming that 
these molecules are not photospheric. Several kinds of 
molecules, with different transition energies, are detected 



in emission. The implication is that these molecules ex- 
tend beyond the background continuum source (the star) . 

In thermal equilibrium, theoretical model atmospheres 
predict that H2O and CO2 are stable in atmospheres 
cooler than 2000-2500 K (i.e., late M) and that SiO is 



stable below 3000 K (Tsuji |196J ). The effect ive te mpera- 
ture of R Set is 4750-5250 K (Shcnton et al. |l994p . Thus, 
these three molecules must be located in 'an extended at- 
mosphere', a region cooler than the photosphere of R Set, 
i.e. located above the photosphere. 

For infrared CO emission, non-LTE effects have 



been proposed (summarized in Oudmaijer et al. 1995). 
However, these non-LTE effects have difficulty explain- 
ing all the emission features, since several molecules with 
various transition energy are seen in emission at the same 
time. It is most likely that the extension of the atmosphere 
around the star is responsible for the emission bands. 

R Set occasionally shows TiO absorption in the visible 
spectra at minimum. TiO is present in the photospheres 
of cool stars, and is used for the spectral classification of 
M-type stars. Its presence has been used to assign R Set 
a spectral class as late as M4 occasion ally a t minimum 
phase (e.g. Preston et al. 1963, Cardelli 1989). The phase 
of the ISO/SWS observation is 0.60 (i.e. intermediate be- 
tween deep minima). It is unlikely that the spectral type 
was early M at the moment of the ISO/SWS observation, 
and there is no evidence that the star ever shows a spec- 
tral type of late M required for stable presence of water- 
vapour. This supports the conclusion that these molecules 
are not photospheric. It seems unlikely that stellar pulsa- 
tion on its own can reduce the stellar temperature from 
- 5000 K to the 3500 K required for the formation of TiO. 
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Fig. 3. (A) The ISO/SWS spectra of R Set (top) are compared with synthesized spectra, which are offset. One 
synthesized spectra includes molecules only; the other also includes a dust excess. Most part of the ISO/SWS spectra 
are well reproduced by the synthesized spectra containing only the molecular bands. A weak dust excess is observed 
at 10 j«m. (B) Synthesized spectra containing only H2O. The black line shows the 2200 K H2O layer only, and the 
grey line shows both H2O layers (2200 K and 1200 K). For comparison, the background source, a 5000 K black body 
is also presented. H2O bands are seen in emission; the global shape of the spectra are dominated by H2O bands with 
an excitation temperature of 2200 K. The 1200 K component contributes to the absorption at 2.7 /im and the emission 
at 6.2 /im. 



Even this molecule may be located in the extended atmo- 
sphere. 

The spectral type of R Set, which is later than that 
of other RV Tau stars, may reflect in part its extended 
atmosphere. The spectral energy distribution is enhanced 
in the infrared region, due to the presence of the non- 
photospheric molecular emission bands in infrared. The 
resultant effective temperature would underestimate the 
actual temperature of the photosphere. Without the ex- 
tended atmosphere, the effective temperature of R Set 
could be more similar to that of other RV Tau stars. 

4.2. The envelope structure of R Set 

Although the molecules are not photosphcric, the excita- 
tion temperatures (800-4000K) and radii (l-10i?*) show 
that these molecules are still located near the star. This 
suggests that they formed recently. 

In our SWS data, the absorption features of H2O 
are seen superimposed on the H2O emission features. 
The cooler layer causing the absorption features and the 
warmer layer causing the emission are therefore seen along 



the same line on sight, with the cooler component mask- 
ing the spectra from the warmer layer. This is possible 
for a disk structure if it is seen approximately edge-on, 
but in general such a configuration is easier to explain if 
the atmosphere consists of an extended shell with a more 
spherical shape. 

The extended atmosphere proposed here describes the 
region above but still close to the photosphere, located 
inside the dust formation radius, at a temperature of 
T ~ 1000 K. The region outside the dust formation ra- 
dius is referred to as the circumstellar envelope. Extended 
molecular atmospheres are seen in pulsating Mira vari- 
ables, where they are caused by the pulsations (e.g., Scholz 
. The extended atmosphere bridges the region be- 
tween the photosphere and the dust forming region, and 
supplies gas from the photosphere to the dust forming re- 
gion. The extended atmosphere and the circumstellar en- 
velope can be distinguished by the temperatures and the 
velocity profiles of the molecules. The molecules in the 
extended atmosphere show complicated velocity profiles, 
often varying with a time scale of one period, and some- 
times moving inward. The circumstellar envelope traces 
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Fig. 4. The ISO/SWS spectra of R Set compared with 
the synthesized spectra (the left scale). Each molecular 
component is normalized to the background source and is 
indicated on the right scale. The quality of the ISO/SWS 
data is lower above 13 fim and the line intensity of each 
CO2 line is not reliable. 



the dust-driven wind, where the velocity can be approxi- 
mated with a near-constant expansion velocity. 

There have been previous indications of this extended 
atmosphere in R Set. Mozurkewich et al. ( 1987 ) found ev- 
idence for velocity variations with a period of 142 days 
in the molecular absorption lines. The variation is at- 
tributed to the pulsation shocks in the extended atmo- 
sphere. This variation is also found in sodium line (Lebre 
& Gillet |1991| ). Hinkle et al. fll996| ) reported that CO 
showed complicated velocity profiles, while OH and H2O, 
which have lower dissociation energies, showed single ve- 
locity profiles. All molecular lines showed time variability. 
They concluded that CO is located near the photosphere 
and is strongly affected by the stellar pulsation, while H2O 
is present in a more extended atmosphere. 

The pulsation shocks propagating into the inner re- 
gion of the molecular atmosphere should not be so strong 
that the molecules are completely destroyed. However, the 
shocks can initially aid the formation of some molecules 



(Duari et al. 1999 ). The molecules may subsequently ac- 
cumulate in a region above the sodium layer. 

We used multiple shells as an approximation to the 
structure of the atmosphere. In reality, a single shell with 
a temperature and density gradient is more likely, where 
the inner, warmer region contributes to the emission and 
the outer, cooler region causes the absorption on the emis- 
sion bands. The temperature and the density will show 
discontinuities at the pulsation shocks, which could still 
lead to multiple discrete shells. However, the shocks do 
not reach 10 and the outer shell is expected to be 
more continuous. 

Recently, H2O molecules have been found in K- and 
M-giants and supergiants whose effective temperatures 



are too high for the pr esence of water (Tsuji |1997| , [2001 
Jennings & Sada 1998| ). These H2O molecules are located 



above atmospheres, as we have found for R Set. The extra 
atmosphere with molecules is commonly seen in cool red 
(super-)giants. 

Many RV Tau stars show an under-abundance of re- 
fractory elements in the photospheric spectra, similar to 
depletion patterns seen in the interstellar medium. To ex- 
plain this effect, re-accretion of gas from a long-lived bi- 
nary disk has been proposed where gas and dust have 



become separated (Waters et al. 1992). Such a disk has 
been reported for the RV Tau star AC Her (Van Winckel 
et al. 1998; Jura et al. 2000p . These disks are located in 
the circumstellar envelope. R Se t shows a relatively small 
depletion (Giridhar et al. 2000 ). Giridhar et al. argued 
that RV Tau stars with late spectral type, such as R Set, 
have enlarged atmospheres: the accreted, depleted gas 
mixes with stellar, non-depleted gas in the extended at- 
mosphere and reduces the depletion. For R Set, there is 
strong evidence that the enlarged atmosphere exists from 
the ISO/SWS observations. 

4.3. Variability 



In Table [A.l| , we summarize published infrared photomet- 
ric data of R Set. One set of observations, taken near 
visible minimum (0=0.07, JD=24433 051), shows fainter 
magnitudes at all infrared bands. Except for this event, 
the magnitudes are comparable at all phases. The magni- 
tudes estimated from the ISO/SWS spectra are also con- 
sistent with other observations at other phases. At the N- 
band, H2O emission increases the flux by a factor of two 
compared to the brightness of the central star (Fig. |J). 
The constancy of the infrared magnitudes suggests that 
the molecules were continuously present around the star 
over the 30 years of observation, and did not form shortly 
before the time of the ISO/SWS observations. 

The visual light curve of R Set shows a number of 
episodes of deep minima over this time, one coinciding 
with the ISO observations and another with the IRAS 
observations (Fig. ^). Before this time, the activity was 
generally less. The extended molecular layers may be re- 
lated to this increased activity, but no spectroscopic or 
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infrared photometric data is available covering the earlier, 
less active period. With an expansion velocity of 5kms _1 , 
the hot water layer lags behind the pulsation of the star 
from the photosphere by less than a year. It is therefore of 
interest that R Set has been relatively quiescent over the 
past 3 years: the effect of this on the water layer could be 
investigated. 



4.4. Mass-loss history 

R Set shows an excess in the IRAS 60 /im band (Fig. |J). 
There is also a weak 12 excess in the IRAS band, 
but our results show that this is mostly due to the H2O 
emission in the 10 /im region, with only a minor con- 
tribution from dust (Fig. |J). The dust mass-loss rate 
which we estimate from the 10 /*m silicate emission is 
Md = 1.5 x 10~ n Mq yr -1 . This mass-loss rate cannot 
produce the observed 60 jjxci excess, and an additional 
cold dust component is needed. 



ISO / sws 

Mode! (molecules + dust) 
IRAS 

Goldsmith etal. (1987) 
4200 K BB (Goldsmith et al.)\ 



Alcolea & Bujarrabal (1991) adopt a two-component 
dust shell model based on the 0.4-100 /mi spectra. For 
the inner shell they derive a dust mass-loss rate of Md = 
4.2 x 10~ n Mq yr -1 , which is comparable to our result 
from the 10 /im excess. For the outer shell they derive 
M d = 2.4 x lO _9 M0yr~ 1 . The CO-derived gas mass- 
loss rate from 12 CO J=l-0 and 2-1 lines is M = 2 x 
10 _7 M Q yr _1 (Bujarrabal et al. |1988| ). If a gas-to-dust 
ratio of 100 is assumed, the gas mass-loss rate derived 
from CO is comparable to the mass-loss rate of the outer 
shell estimated from the dust shell model. Also, the low 
excitation temperature of the millimeter CO rotation line 
suggests that the CO traces the past mass-loss rate. The 
implication is that the mass-loss rate of R Set has reduced 



by about a factor of 100 (Alcolea & Bujarrabal 1991). 
Using the inner radius of the outer (cold) dust shell of 
3.7 x 10 16 cm and the CO expansion velocity of 5 kms" 1 , 
Alcolea & Bujarrabal estimated the elapsed time after the 
large mass loss ceased as ~ 2000 yr. This issue will be 
discussed in the next section. 

The low current dust mass-loss rate raises a peculiar 
point. The SiO column density derived here is compara- 
ble to that of Mira variables (Yamamura et al. 1999| ). 
R Set should therefore have a sufficiently large amount 
of SiO for silicate dust formation. However, the current 
dust mass-loss rate is two orders of magnitudes lower than 
that of Mira variables (typically M<j = 10~ 9 Mq yr _1 - 
10~ 8 Mq yr -1 ). This suggest a process to inhibit the dust 
formation from the molecules. The hotter radiation field 
from the star, compared to Miras, may affect the dust 
formation sequence. 

For AGB stars, the extended atmosphere is thought 
to be built up by pulsations, while radiation pressure on 
dust is the main mechanism driving the mass loss (e.g., 
Fleischer et al. |1992| Hofner & Dorfi |1997| ). It is there- 
fore conceivable that both low metallicity and less dust 
grains operate at small mass-loss in R Set, even if the in- 




1 10 
Wavelength [|j,m] 



Fig. 5. The ISO/SWS spectra are shown with the IRAS 
fluxes (the 100 /jm band is an upper limit) and UV and IR 
photometric data (Goldsmith et al. |1987| ). A 4200 K black 
body, fitting the UV and IR photometric data (Goldsmith 
et al. 1987) is also plotted. Zero points and effective wave- 



length for the filters from Evans (1994) are adopted for 
the data in Goldsmith et al. (1987) 



termediate reservoir of SiO has the same size as in Mira 
variables. 

4.5. Evolutionary status 

RV Tau stars are commonly assumed to be post-AGB 
stars with low initial mass. For the RV Tau stars with 
IRAS 60 /im excess, the cessation of the mass loss oc- 
curred about five hundreds years ago (Jura 1986 ). R Set 
also shows evidence for a larger mass-loss rate in the past, 
and is also sug geste d to be in the post-AGB phase (Alcolea 
& Bujarrabal |l99i| ). 

The mass loss ceased at the time the star would have 
left the AGB. The effective temperature (T e e) at AGB 
is in the range of 2000-3000 K, and the current effec- 
tive temperature of R Set is in the range 4000-6000 K. 
Using the time since the mass-loss reduction of 2000 yr, 
we obtain a rate of increase of the stellar temperature 
dT c ff/dt = 1-1.5 Kyr - 1 . Th is rate can be compared to 
calculations by Blocker ( |l995|) of the temperature gradient 
as function of the post-AGB evolutionary time scale (in 
terms of log di/dT G ff yr K _1 , Fig. 6 in his paper). The post- 
AGB time is calculated from the end of the high mass-loss 
rate. The rate we find for R Set is reproduced for the early 
post-AGB evolution by the model tracks of stellar (core) 
mass range of 0.61, 0.57 and 0.55 Mq. These masses are at 
the lower end of post-AGB mass ranges, and may indicate 
a progenitor star with low initial mass. 

Together with the evolution towards higher T e g, the 
period of the star will also reduce as the stellar radius 
becomes smaller. R Set is one of the oldest variable stars 
known, discovered in 1785 by Edward Pigott, observing 
from York. He observed the star regularly until 1803. Later 
it was observed extensively by Argelander from Bonn, be- 
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tween 1843 and 1862, and very good coverage is available 
from this time. This long-term light curve allows us to 
search for evidence of period evolution. 

If all RV Tau stars are post-AGB, they are expected 
to evolve at constant luminosity towards higher effective 
temperatures. We use the period equation for AGB stars 



logP = L9491ogii-0.91ogM-2.07 



(1) 



from Wood (1990), where P is the period in days, R and 
M are the stellar radius and mass in solar units. Taking 
L = 4ttR 2 T^ s = const, where L is the luminosity, we can 
rewrite the equation as 



logP = -3.898 log T eS + 0.975 log L- 0.9 log M- 2.62. (2) 

This implies that the rate of period change is closely 
related to the temperature evolution: 



P 
P 



-3.9- 



off 



(3) 



The right-hand side corresponds to the calculated tem- 
perature increase of ~ 1 Kyr -1 . Using the present period 
of 140 days, a change of —11 days per century is expected. 
This should be easily detectable in t he am ateur archives 
and historical observations. Zsoldos ( 1995 ) has noted the 
lack of observed secular period changes in other RV Tau 
stars and stressed that the low-mass post-AGB interpre- 
tation must be considered uncertain. 

We have analyzed published observations of R Set 
to search for evidence of period decrease. Individual vi- 
sual observations were taken from the databases of the 
AFOEVj], the BAAVSSg and the VSOLjf], which extend 
back over more than a century. The top panel of Fig. ^| 
shows the light curve after averaging into 5-day bins. The 
deep minima occur with a period of about 140 days, and 
they alternate with the much shallower second ary m in- 
ima. We used wavelet analysis (Bedding et al. 1998 ) to 
measure period evolution in both the 140-day period and 
its 70-day harmonic. These are shown in the middle and 
bottom panels of Fig. [| respectively. While R Set shows 
substantial period jitter, there is no evidence for a secular 
trend in the mean period during the past century. 

The evolution of the primary and harmonic periods is 
plotted in Fig. |] as open symbols. We also show closed 
symbols which are derived from wavelet analysis of earlier 
published dates of minima. Those data are sparse, with 
some measurements around 1800 and then a large gap 
until 1860. 

Period determinations during the 1800s were also made 
by the observers themselves. The values are somewhat 
scattered, but fall into two groups: around 60 days and 
around 70 days. A summary of t hese d eterminations can 
be found in Muller & Hartwig ( 1918 ). As they discuss, 



Pigott found periods of 63 and 61.5 days, based on the pre- 
1800 observations. Westphal found 68.8 and 60.56 days, 
both determinations based on the same data! Argelander 
found 60 days in 1844, based on one year of observa- 
tions, but later derived 71 days from a 15-year light curve. 
Schmidt, in 1879, was the first to derive the fundamental 
period and found 140 days. 

All determinations before Schmidt clearly refer to the 
harmonic period. The reason is that the observers mea- 
sured the times of maxima and minima and used these 
in their analysis: the primary and secondary minima were 
not distinguished in the published observations. Since the 
two corresponding frequency groups differ by about 1 cy- 
cle per year, it seems probably that one of these periods 
is an alias of the other. The current (harmonic) period 
is about 70 days (e.g., 70.6 days measured by Pollard et 
al. 1996 ). It therefore seems that the (harmonic) period 
has been about 70 days for as long as it has been mea- 
sured. (Ascribing reality to the 60-day period would fur- 
ther worsen the discrepancy with eq. (^), since it would 
give a period evolution with the wrong sign.) 

We re-reduced Argelanders observations (1869) to 



check his period determination. The Fourier power spec- 
trum shows a significant peak at 145.0 days, plus numer- 
ous weaker peaks in two ranges, centred at 140 and 70 
days. Thus we confirm his period, but it is clear that the 
true dominant period was already evident in these data. 
This analysis also confirms the period stability of R Set. 

In conclusion, over the past two hundred years there 
is no evidence for a change of more than 2 days in the 
harmonic period, corresponding to 4 days in the main pe- 
riod. This limit is a factor of 5 below the predictions from 
Eq. (§). 



Percy et al. (1991) analyzed 150 years of data for R Set 
and reported that a sudden decrease of 1.0 d occurred in 
the harmonic period in 1941. However, their calculation of 
0.07 day/year as the implied average rate of period change 
during the 150 years of observation is incorrect. The cor- 
rect value is ten times smaller (1/150 = 0.007day/year). 
Thus, their conclusion that R Set shows a period decrease 
in good agreement with evolutionary models is not valid. 

The lack of evidence for period evolution implies that 
R Set is presently not evolving towards higher tempera- 
ture. This is in contrast to the post-AGB tracks and sug- 
gests alternative evolutionary phases should be explored, 
in agreement with the discussion by Zsoldos ( 1995 ). 

The past mass-loss rate of R Set is of order 
10~ 7 MQyr" 1 . This is rather low for the final mass-loss 
rate of the AGB, which is expected to be more than 
~ 10" 5 M yr" 1 (Bowen & Willson [l99l] ). In addition, the 



abundance of s-process elements for R Set is low. One pos- 
sible explanation is that R Set may not yet have reached 



1 http:/ /cdsweb.u-strasbg.fr/afoev/english.htx 

2 http://www.britastro.org/vss 

3 http:/ /www. kusastro.kyoto-u.ac.jp/vsnet/VSOLJ/vsolj. html 



the tip of the AGB (Giridhar et al. |2000J). 

It is therefore possible that R Set is still an AGB 
star. A sudden decrease of the mass- loss rate on the AGB 
is expected following the thermal pulse: the lower lumi- 
nosity during the quiescent helium burning gives both a 
higher temperature and a lower mass-loss rate (Zijlstra et 
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al. |1992| , Vassiliadis & Wood |1993| , Blocker ]! 995Q . The pe- 
riod evolution slows down and reverses during the helium- 
burning phase. The phase lasts for ~ 5000 yr. This phase 
may fit the period evolution, and time scale for the de- 
tached dust shell, quite well. 

We therefore suggest the possibility that R Set is not a 
post-AGB star, but is an AGB star in the helium-burning 
phase of the thermal-pulse cycle. 

5. Conclusions 

We find H 2 0, SiO, C0 2 , and CO in the ISO/SWS spec- 
tra of R Set. H2O emission dominates the near- and mid- 
infrared spectra. The molecules are probably located in a 
spherical extended atmosphere. Pulsations may gradually 
lift up the gas from the photosphere to the extended at- 
mosphere. Molecular formation appears to take place in 
the extended atmosphere. 

The column density estimated from the infrared SiO 
bands is comparable to that in Mira variables, however, 
the current dust mass-loss rate of R Set is lower than Mira 



variables by two orders of magnitudes. The high effective 
temperature may prevent dust formation around this star. 

While a high mass-loss rate in the past is indicated by 
an excess at 60 /Ltm, and by the CO millimeter wavelength 
emission, the current mass-loss rate is two orders of magni- 
tude lower: the dust mass-loss rate M<j = 10~ n Mq yr -1 . 
The reduction in the mass-loss rate, and the time scale, 
can be explained if R Set is a low-mass post-AGB star. 
However, the period evolution is too slow for post-AGB 
evolution. We therefore suggest that R Set may be an AGB 
star in the helium-burning phase of the thermal-pulse cy- 
cle. 

As future work, monitoring of the infrared photomet- 
ric bands could help in understanding the structure of the 
extended atmosphere in this star, as well as in Mira vari- 
ables. The infrared light curve of Mira variables shows a 
phase delay of maximum phase by about </>=0.1. For the 
infrared variability of Mira variables, dust might be im- 
portant (e.g. Winters et al. 1994), but molecules in the 



pulsating atmosphere contribute as well (Matsuura et al. 
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( |1999D . The zero point for the 10.4 /iin filter is calculated 
from the spectrum of Vega (Cohen et al. 1995). 
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Fig. 7. The period variation of R Set over 200 years. 
The upper figure is for the primary period and the lower 
figure for the harmonic period. The upper horizontal scale 
shows the years. The predicted period variation from the 
post- AG B evolutionary track is shown by the dotted lines 
(Sec. 15). The observed variation is shown by the dashed 
line, where the harmonic period is fitted with the linear 
function of P = [70.71 ± 0.33] - [0.03 ± 0.10] x 10~ 4 x 
(JD — 240 000): i.e., the period is constant over 200 years 
(long dash line). 



2002j ) . If the same infrared phase delay is found for R Set 
the contribution of molecules on t he infrared light-curve 
will be strengthened. In Table A.l, the infrared photom- 
etry shows slightly fainter magnitudes at phase <p =0.07 
(JD=244,3051.31). More frequent observations of the in- 
frared light curve could find a phase delay. 
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Appendix A: Infrared photometry 

Infrared photometric data taken from the literature are 
The data are taken from Gehrz ( 1972j ), 



listed in Table A.l 



Lloyd Evans ( |1985| ), Goldsmith et al. ( |l987j ), and Shenton 
et al. (1994). The zero points for the IRAS bands are 
taken from Wainscoat et al. ( 1992] ). The flux density of 
the ISO/SWS spectra is converted to magnitude, using 
filter transmission in van der Bliek et al. (1996) and zero 
points in Evans (1994) the for L- and M-bands. For N- and 
Q-bands, effective centre wavelengths and widths in Evans 



(1994) are used. We measured an N-band magnitude with 
TIMMI-2 on 3.6 m Telescope of ESO on 4 August, 2001, 
using a narrow band filter at 10.4 /im. HD 133774 was used 
as a calibration source adopting the template in Cohen 
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Table A.l. Infrared photometric data 



JD 
2440000+ 


Phase 


J 


H 


K 


L 


M 


N 


Q 


IRAS 12 IRAS 25 


Ref 


Uncertain* 








2.8 


1.8 


1.4 


0.9* 






1 


2682.31 


0.64 


2.88 


2.33 


2.04 


1.53 










2 


2999.36 


0.72 


2.78 


2.28 


2.02 


1.51 










2 


3051.31 


0.07 


3.38 


3.04 


2.80 


2.33 










2 


6243.61 


0.61 


2.84 


2.26 


2.00 


1.49 










3 


6248.53 


0.64 










1.50 


0.85 






3 


6253.50 


0.68 


2.84 


2.26 


2.00 


1.49 










3 


6991.48 


0.80 


3.10 


2.52 


2.26 


1.77 


1.15 


0.36 


0.46 




4 


7000.49 


0.86 


3.10 


2.49 


2.22 


1.65 










4 


7040.25 


0.15 


3.11 


2.57 


2.36 


1.91 


1.67 


0.93 


-0.01 




4 


8840.36 


0.49 


3.06 


2.46 


2.26 


1.91 










4 


8895.36 


0.86 


2.84 


2.31 


2.10 


1.76 










4 


Uncertain 


















0.74 s 0.03 § 


5 


10153.24 


0.60 








1.81* 


1.36 1 


0.98 1 


0.00 1 


0.75 s 0.05 s 


6 


12125.55 














0.75 
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Shenton et al. dl994[ ). 5: IRAS data. IRAS observed R Set 
Es timat e from the ISO/SWS data. Filter transmission of 



1: Gehrz ( |1972[ ). 2: Lloyd Evans fll985| ). 3: Goldsmith et al. ( |1987| ). 4 
three times, and no variabil ity h as been found (Var flag is %). 6: 
van der Bliek et al. (L, M; 1996) and IRAS Ex planat ory Supplement (1988) are adopted. N and Q band are estimated from 
effective centre wavelength and width in Evans (1994). 7: ESO/TIMMI-2 observations. Not broad JV-broad band but narrow 
TV-band filter (#10.4). 

^ Observed b etwee n April 1970 and April 1971. * Using a narrow band filter centred at 10.8 /jm. * Zero magnitudes are taken 

from Evans (1994). s Zero magnitudes are taken from Wainscoat et al. (1992). 

Phases from JD=2442682.31 to JD=2446253.50 are estimated from the light curve in Mattei ([t985|; |l99(]). 
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